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^ ' In the standard hot big bang theory, when the Universe was about 1 — 10 /is old, the cosmological 

matter is conjectured to undergo Quantum Chromodynamics (QCD) phase transition(s) from quark 
matter to hadrons. In the present work, we study the cosmological quark-hadron phase transition 
in two different physical scenarios. First, by assuming that the phase transition would be described 
by an effective nucleation theory (prompt first-order phase transition), we analyze the evolution 
of the relevant cosmological parameters of the early Universe (energy density p, temperature T, 
Hubble parameter H and the scale factor a) before, during and after the phase transition. To study 
the cosmological dynamics and the time evolution, we use both analytical and numerical methods. 
The case where the Universe evolved through a mixed phase with a small initial supercooling and 
monotonically growing hadronic bubbles is also considered in detail. The numerical estimation of the 
cosmological parameters, a and H for instance, shows that the time evolution of the Universe varies 
^ ^ from phase to phase. As the QCD era turns to be fairly accessible in the high-energy experiments 

and the lattice QCD simulations, the QCD equation of state is very well defined. In light of these 
QCD results, we develop a systematic study of the crossover quark-hadron phase transition and an 
estimation for the time evolution of the Hubble parameter during the crossover . 

PACS numbers; 04.50.-h, 98.80.Cq, 98. 80. Bp, 98.80.Jk 



I. INTRODUCTION 

According to the standard model of cosmology, as the Universe extremely expanded and cooled down, it is 
likely to expect that the cosmological background matter should undergo a series of symmetry-breaking phase 
transitions, at which various topological defects may have formed. The study of phase transition from quark- 
gluon plasma (QGP) to hadrons in the early Universe dates back to about three decades ago [ln5]. A first-order 
phase transition in various scenarios is assumed to take place @. In one scenario, it has been suggested that QGP 
thermodynamically condensates into a hadron gas. In the second scenario, it is conjectured that the Universe 
was being supercooled and an out-of-equilibrium nucleation of hadron bubbles in the QGP surrounding should 
take place. In the third scenario, it has been argued that the phase transition took place in accompany with 
a small supercooling. Apparently, the coexistence of hadrons and QGP is accessible after the nucleation. The 
latter would generate fluctuations in the isothermal baryon density i.e., inhomogeneity and therefore can lead to 
drastic astrophysical consequences. From Yang-Mills theory, we have learned a lot about the kinetics and order 
of the phase transition The lattice QCD is a reliable method describing the strongly interacting matter 
for the whole temperature range starting from very low temperatures (ground state) to very high temperatures 
(perturbative QCD). Recently, a remarkable discovery of the QGP properties has been achieved in the heavy-ion 
collision program [8l-[ll|. The QGP is likely a strongly correlated phase with finite bulk and shear viscosity. 

A first-order phase transition is proceeded by bubble nucleation and rapid expansion. When at least 4 — n of 
these bubble collide, where n = 0, 1, 2, an n-dimensional topological defect may form in the region between them 
. Recent lattice QCD calculations for two quark flavors suggest that QCD reliably describes a transition at 
Tc ^ 173 McV [l^. It is neither first- nor second-order. With increasing temperature there is a rapid change 
in all thermodynamic quantities. This phase transition, which could have occurred in the early universe, could 
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lead to the formation of relic quark-giuoii plasma objects, which still survive today. It will be elaborated below, 
that the order of the phase transition strongly depends on the mass and flavor of the quarks. 

As given above, studying the first-orde r q uark-hadron phase transition in the early Universe has a long 
history. It can be characterized as follows [1^. As the color deconfined QGP cools down below Tc, it becomes 
energetically favorable to form color confined hadrons (primarily the lightest Goldstone bosons; the pions and 
a tiny amount of neutrons and protons, due to the conserved net baryon number). However, the new phase 
does not show up immediately. A characteristic feature of the first-order phase transition is that a part of the 
supercooling is needed to overcome the energy expense of forming the surface of the bubble and the new hadron 
phase. When a hadron bubble is nucleated, latent heat is released and a spherical shock wave expands into the 
surrounding supercooled QGP. This reheats the plasma to the critical temperature, preventing further nuclcation 
in a region passed by one or more shock fronts. Generally, the bubble growth is described by deflagrations with 
a shock front preceding the actual transition front. The nuclcation stops, when the whole Universe has reheated 
to Tc- This part of the phase transition passes very fast, in about 0.05 //sec, during which the cosmic expansion 
is totally negligible. After that, the hadron bubbles grow at the expense of the quark phase and eventually 
percolate or coalesce. When neglecting the possibility of the quark nugget production, the transition is assumed 
to stop, when all QGP has been converted to hadrons. 

Depending on the numerical values of the parameters, both deflagrations and detonations can appear. The 
hadron bubbles can nucleate at very large distance scales and the phase transition may be completed without 
reheating to the critical temperature. During the low temperature phase in the phase transition the bubble 
can grow as a supersonic deflagration consisting of a Jouguet deflagration followed bil^ rarefaction wave. The 
velocity of the supersonic deflagration varies between the sound and light velocities [IJ] . The small-scale effects 
of flnite wall width and surface tension have been incorporated in a numerical code, also including both the 
complete hydrodynamics of the problem and a phenomenological model for the microscopic entropy production 
mechanism at the phase transition surface [l^ . The decaying droplets leave behind no rarefaction wave, so that 
any baryon number inhomogeneity generated previously should survive the decay. 

The nuclcation of bubbles, the collisions of shock fronts preceding the bubble, the arrestation of the bubble 
growth by the reheating, the condensation of the baryon number and the resulting density perturbations af- 
ter a first-order phase transition through the mixed phase have been studied in a scenario with small initial 
supercooling and monotonically growing hadronic bubbles [Hj . The growth of bubbles after the initial nucle- 
ation event in the generic first-order cosmological phase transitions, which is characterized by the latent heat 
L, the interface tension a and the correlation length ^ and is driven by a scalar order parameter (f> has been 
considered in Ref. [l^. The mean distance of the nuclcation d^uc in a first-order cosmological quark-hadron 
phase transition has been introduced in Ref. [l7| . For a homogeneous nuclcation dnuc ^ 2cm. On the other 
hand, the impurities can lead to heterogeneous nuclcation, with dnuc of several meters. The latter value could 
change the outcome of the big bang nucleosynthesis. The study of the hydrodynamics of the disconnected quark 
regions during the final stages of the cosmological quark-hadron transition has been carried out in Ref. [l^ . 
It has been shown that a self-similar solution likely exists. The inclusion of the relativistic radiative transfer 
produces significantly different results. Furthermore, it enables the formation of high density regions at the end 
of the drop evaporation p^ . The linear stability analysis of the relativistic detonation fronts, representing the 
phase interface in first-order phase transitions, has been performed in Ref. [23|. The strong detonations are 
evolutionary and stable with respect to the corrugations of the front. Moreover, Chapman- Jouguet detonations 
appear to be unconditionally linearly stable. Taking into account the simultaneous effects of the baryon number 
flux suppression at the phase interface, the entropy extraction by means of the particles having long mean free 
paths and baryon diffusion shows that signiflcant baryon number concentrations, up to densities above that of 
nuclear matter, represent an inevitable outcome within this scenario [20j . 

The abundance and size distribution of the quark nuggets formed a few microseconds after the big bang due 
to a first-order QGD phase transition have been estimated in Ref. [2l|. The evolution and the collision of 
slow-moving true vacuum bubbles are examined in Ref. [2^ . The comoving bubble walls prevent the formation 
of extra defects and may lead to an increase of any primordial magnetic field. Within an effective model of 
QGD, the quark-hadron phase transition was studied in Ref. [1^. In a reasonable range of the parameters 
of the model, bodies with a quark content between 10^^ and 10 Mq could have been formed in the early 
universe. A significant amount of entropy is released during the transition. The density fluctuations amplified 
by the vanishing sound velocity effect during the quark-hadron phase transition could lead to QGP lumps 
decoupled from the expansion, which rapidly transform to quark nuggets (2^ 1 . The inhomogeneous nuclcation, 
as a new mechanism for the cosmological QGD phase transition, was proposed by Ignatius and Schwartz [2^. 
In this model the typical distance between bubble centers is of the order of a few meters. The resulting baryon 
inhomogencities may affect the primordial nucleosynthesis. 



3 



Recent lattice QCD simulations turn to be able to provide an accurate tool to study - among others - the 
thermodynamics of the strongly interacting matter. The critical temperature was a subject of different lattice 
QCD simulations [2^-[13- know so far that for two quark flavors {nj = 2) the transition is second-order or 
rapid crossover and Tc ~ 173±8 MeV. For n/ = 3, we have a first-order phase transition and Tc — 154±8 McV. 
For 71/ = 2 + 1 i.e., two degenerate light quarks and one heavy strange quark, the transition is again crossover 
and Tc ~ 173 ± 8 MeV. For the pure gauge theory, Tc — 271 ± 2 MeV and the deconfinement phase transition is 
first-order. In all these lattice QCD simulations, the quark masses are much heavier than their physical values. 
With recent computational facilities and modern algorithms, it is now possible to use values very close to the 
physical masses. This raised the critical temperature, for instance, Tc ~ 200 MeV for = 2 -f 1. From this 
discussion, wc conclude that the order of the phase transition can be cither continuous or discontinuous. It 
depends - among others - on the quark flavors and their masses. The extreme conditions in the early universe, 
like high temperatures, high densities and out-of-thermal and out-of-chemical equilibrium, likely affect the 
properties of the partonic matter. Yet, we have no access to study this issue. Recent lattice QCD outputs 
have been used in j33j to work out the expansion law of the Universe during the cosmological quark-hadron 
transition. The cosmological behavior found using lattice data was compared with the one obtainable in case 
the transitions were first-order . The differences between these two scenarios are too small to be tested with 
cosmological data, but the coming of the era of precision cosmology might open the possibility of testing the 
nature of the QCD transition by using cosmological data. 

In the present work, we consider two cases. First, we assume that the phase transition is of first-order. The 
cosmological evolutions during the quark and hadron phases are investigated in detail. The main cosmological 
parameters are obtained for each phase. The hadron fraction ft,, whose time evolution describes the conversion 
process, is an important parameter to describe the phase transition and its expression is obtained in an analytical 
form, h seems to behave as an order parameter. The second part of this study is devoted to an extension of 
previous works [34| - |40j . in which we have applied the equations of state deduced from recent lattice QCD 
simulations at almost physical masses and more accurate lattice configurations in order to study the cosmology 
of the early universe. With the use of these equations of state we can study the evolution equations of the main 
physical parameters of the cosmological models. In light of these QCD results, we develop a systematic study of 
the crossover quark-hadron phase transition and an estimation for the time evolution of the Hubble parameter 
during the crossover in the presence of bulk viscous effects. 

This paper is organized in the following manner. In Section|TTl the background geometry and the gravitational 
field equations are written down and the description of the viscous effects in different theoretical models is 
presented. In Section Hill wc lay down the equations of state and the relevant physical quantities, necessary for 
the discussion of the first-order quark-hadron phase transition. In Section HVl wc analyze in detail the dynamics 
of the Universe during first-order quark-hadron phase transition. The phase transition in the lattice QCD 
simulations and the heavy-ion collisions and the QCD equation of state (EoS) are discussed in Section |Vl The 
cosmological evolution of the Universe during the crossover in the presence of bulk viscous effects is analyzed in 
Section The cosmological implications of our results are discussed in Section IVIII We discuss and summarize 
our results in Section IVIIII 

In the present paper we use natural units with c=/i = fcB = 1, in which 87rG = 1/mp; = 1.687 x 
10-43 MeV"^, where mpi is the "reduced" Planck mass. The "reduced" Planck time is given by ip; = l/nipi = 



We assume that the early Universe is filled with a bulk viscous cosmological fluid and its geometry is given 
by a spatially fiat Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 



where a{t) is the dimensionless scale factor, which describes the expansion of the universe. At a vanishing 
cosmological constant, the Einstein gravitational field equations in the flat Universe read 



4.0 X 1152 X 10-22 MeV-i. 



II. GEOMETRY AND FIELD EQUATIONS 
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The energy-momentum tensor of the bulk viscous cosmological fluid filling the very early Universe is given by 

m 

^{p + p + n)u,u'' ~ip + Il)S'[, (3) 

where indices i,k take discrete values 0,f,2,3, p is the energy density, p is the thermodynamic pressure, 11 
is the bulk viscous pressure and Ui is the four velocity, satisfying the normalization condition UiU^ = 1. The 
particle and entropy fluxes are defined according to iV* = nu* and S*' = sA^* — (tIP /2£^T^ w*, where n is the 
number density, s is the specific entropy, T is the finite temperature, ^ is the bulk viscosity coefficient and r 
gives the relaxation coefficient for the transient bulk viscous effect (i.e. the relaxation time), respectively. The 
evolution of the cosmological fluid is subject to obeying the dynamical laws of the particle number conservation 
N'\^ = and Gibbs' equation Tdp ~ d{p/n) +pd{l/n) [4f|. In the following, we shall also suppose that the 
energy- momentum tensor of the cosmological fluid is conserved, i.e., T^f, = 0, where ; denotes the covariant 
derivative with respect to the metric. 

The bulk viscous effects can generally be described by means of an effective pressure 11, formally included in 
the effective thermodynamic pressure Pef f = p + H [4l| . Then in the comoving frame the energy-momentum 
tensor has the components Tq = p, = r| = Tg = —Pef / ■ For the line element given by Eq. ([1]) , the Einstein 
field equations read 

H + H' = --Xr (ip.ii + p). (5) 

where one dot denotes the derivative with respect to the time G is the gravitational constant and H{t) = 
d(t)/a{t) is the Hubble parameter. Expressions Q and ([5]) lead to a generic expression for the time evolution 
of iJ: 

H = --V(p + p + n). (6) 

From the field equations or with the use of the conservation of the energy-momentum tensor we obtain the 
following equation (Bianchi identity), relating the time variation of the energy density to the Hubble parameter: 

p + 3{p + Peff)H^0. (7) 

In order to solve the field equations, we necessarily need an equation of state and an estimation for the bulk 
viscous H, characterizing the viscous properties of the matter in the expanding universe. 



A. Eckart relativistic viscous fluid 

The first attempts at creating a theory of relativistic fluids were those of Eckart [i^] and Landau and Lifshitz 
(43I . These theories are now known to be pathological in several respects. Regardless of the choice of the equation 
of state, all equilibrium states in these theories are unstable and in addition signals may be propagated through 
the fluid at velocities exceeding the speed of light c violating the causality principle. These problems arise due 
to the nature of the first-order of this theory [review Eq. ©], that it considers only the first-order deviations 
from the equilibrium leading to parabolic differential equations, because of the infinite speeds of propagation for 
the heat flow and viscosity, which contradicts the principle of causality. Conventional theory is thus applicable 
only to phenomena which are quasistationary, i.e., slowly varying on space- and time-scales characterized by 
mean free path and mean collision time. 

The Eckart theory can be applied on modelling the cosmic background fluid as a continuum with a well- 
defined average 4- velocity field it" where u^Ua ~ —1. The vector number density = nu°' can be estimated, 
when unbalanced creation/annihilation processes take place; n"^ = 0. This apparently means that 

n + 3Hn = 0, (8) 

where the Hubble parameter H = u"^. In the case of a viscous fluid, the entropy current. 



b ~ snu , 



(9) 
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is no longer conserved. The covariant form of second law of thermodynamics is S"^ > and the divergence of 
entropy current is given by TS°^ = —SHU. This is another feature of Eckart's theory. It violates the second 
law of thermodynamics. 

The evolution of the cosmological fluid is subject to the dynamical laws of particle number conservation 
iV.*j = and Gibbs' equation Tdp = d{p/n) +pd{l/n). Then, from the Gibbs equation, the covariant entropy 
current can be obtained as 

U=-3^H. (10) 

This is a linear first-order relationship between the thcrmodynamical flux IT and the corresponding force H. 
Substituting in Eq. ((6)) results in 

H = -^iP + P~^^H). (11) 



B. Israel-Stewart relativistic viscous fluid 

A relativistic second-order theory was introduced by Israel and Stewart [3 EE] and further developed by 
Hiscock and Lindblom j46j through the exienrferf irreversible thermodynamics. In this model, the deviations from 
equilibrium (bulk stress, heat flow and shear stress) are treated as independent dynamical variables, resulting 
in 14 dynamical fluid variables to be determined. The causal thermodynamics and its role in general relativity 
are reviewed in Rcf. [U 
thcrmodynamical flux 11 



A general algebraic form for including a second-order term in the dissipative 
44143 reads 



S^ = snu^+(3n'—, (12) 

where /? is a proportionality constant. 

For the evolution of the bulk viscous pressure, we adopt the causal evolution equation [4l| obtained in 
the simplest way (linear in 11) to satisfy the 7?-theorcm (i.e., for the entropy production to be nonnegative, 
S.^j = iP/^T > m, 14511 '). According to the causal relativistic IS theory, the evolution equation of the bulk 
viscous pressure reads |4l| 

Ttl + n = -3^H-^rIl(^iH+^-^~^y (13) 

where r is the relaxation time. In order to have a closed system from Eqs. (O and ([T^ . we have to take into 
consideration equations of state for the pressure p, the temperature T and the relaxation time r, respectively. 



III. FIRST-ORDER QUARK-HADRON PHASE TRANSITION 



In this Section, we outline the relevant thermodynamic quantities of the quark-hadron phase transition, which 
will be used in the following sections. Note that the scale of the cosmological QCD transition is given by the 
Hubble radius Rh at the transition: Rh ~ mpi/T^ ^ 10 km, where Tc is the critical temperature. The mass 
inside the Hubble volume is ^ 1 AIq. The expansion time scale is 10~^ s, which should be compared with the 
time-scale of QCD, 1 fm/c~ 10"^'^ s. Even the rate of the weak interactions exceeds the Hubble rate by a 
factor of 10^. Therefore, in this phase the photons, the leptons, the quarks and the gluons (or pious) are lightly 
coupled and may be described as a single, adiabatically expanding fluid [T7| . 

At high temperatures T > Tc, the baryon number density ub may be defined as ub = (1/3) ('^g ~ 
where riq (riq) is the number density of a specific quark (anti-quark) flavor and the sum is taken over all quark 
flavors. In utilizing these relations, it is apparent that QGP matter is assumed to be characterized as an ideal 
gas. At r < 1 GeV only the u, d and s quarks contribute significantly. At low temperatures T < Tc the baryon 
number density is defined as = J2 ("^6 ~ ^s)? with the summation extended over all baryon species b. In 
order to study the quark-hadron phase transition it is necessary to specify EoS of the matter, in both quark and 
hadron state. Giving an equation of state is equivalent to give the pressure as a function of the temperature T 
and chemical potential fi. 
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At high temperatures the quark chemical potentials are equal, because of the weak interactions which appar- 
ently keep them in chemical equilibrium and the chemical potentials for leptons arc assumed to vanish. Thus 
the chemical potential for a baryon is defined hy fiB = 'ifJ-q- The baryon number density of an ideal Fermi gas 
of three quark flavors is given by ub — T^y^ts/S, leading to hb/T ~ 10~^ at T > Tc- At low temperatures 
Hb/T ^ 10^^. Therefore the assumption of a vanishing chemical potential at the phase transition temperature 
in both quark and hadron phase represents an excellent approximation for the study of EoS of the cosmological 
matter in the early universe. In addition to the strongly interacting matter we assume that in each phase there 
are present leptons and relativistic photons, satisfying equations of state similar to that of hadronic matter [l^ . 



A. Thermodynamic parameters of the quark and hadronic matter 

The equation of state of the ideal gas in QGP phase can generally be given in the form 

^ 3agT^ + V{T), (14) 
Pg - agT^^V{T), (15) 

where V (T) is the self-interaction potential. Uq = {tt'^/90) gq, with .g, = 16 + (21/2)iV_F + 14.25 51.25 and 
Np = 2. As given in Ref. [2^, the self-interaction potential reads 

V (T) = B + "/tT^ -arT^, (16) 

where B is the bag constant, ax = 77r^/20 and 77- = nil/A, with nis is the mass of the strange quark G (60 — 200) 
MeV. The form of the potential V corresponds to a physical model in which the quark fields are interacting 
with a chiral field formed with the tt meson field and a scalar field. If the temperature effects can be ignored, 
EoS in the quark phase takes the form of the MIT bag model equation of state, Pq = {pq — 4i?)/3, MIT stands 
for Massachusetts Institute for Technology. The results obtained in the low energy hadron spectroscopy, the 
heavy- ion collisions and the phenomenological fits of the light hadron properties give an estimation for i?^/^. It 
ranges between 100 and 200 MeV [iTj . 

In the hadron phase, we assume that the cosmological fluid is consisting of an ideal gas of massless pions 
and nucleons described by the Maxwell-Boltzmann statistics. The energy density p^ and pressure ph can be 
respectively approximated by 

Ph{T) = ^ph{T)^a^T\ (17) 

where = (tt^/QO) gt and gt — 17.25. For the entropy densities s{T) = dp/dT in the two phases we obtain 

s,{T) - -27TT + 4(a,-HaT)T3, (18) 
Sh[T) = Aa^T\ (19) 

The critical temperature Tc is dcflned by the condition pq (Tc) = ph (Tc) [l2l | and is given, in the present model, 

by 



J.2 ^ 7T + Vtt + 4:{aq + aT- a^) B ^^^^ 
^ 2(aq + Q;r — a^) 

For nis = 200 McV and B^/'^ = 200 MeV, the transition temperature is of the order of Tc ~ 125 MeV. 
According to the first-order of the phase transition, all the physical quantities, like the energy density, pressure 
and entropy, exhibit discontinuities across the critical curve. At the critical temperature, the ratios of the 
relevant physical quantities, the energy and the entropy density, respectively, are given by 

Pq {Tc) ^ 4:aqT;^ - Pq (Tc) ^ ia^T^ - ph {Tc) ^ 4aq - 
Ph {Tc) 3a^r4 3a^r4 3a^ 

and 



(Tc) _ TtOtt + (og + ar) ^/^^ + 4:{aq + ar - a^) B 

[ZZ) 



{Tc) (^j^ + ^/-f'^ + 'i{aq + aT-a^)B^ 
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respectively. For = 200 MeV and B^^^ = 200 MeV, the ratios pq (Tc) /ph (Tc), given by Eq. ^ and 
(^c) /sh (Tc), given by Eq. (f^ . equal 3.62 and 4.628, respectively. So far, we conclude that the energy 
density and entropy suddenly decrease to nearly one-fifth of its value, when the system undergoes a first-order 
phase transition at Tc- According to the first-law of thermodynamics, the entropy s can be expressed in terms 
of the pressure p and the energy density p, so that at vanishing chemical potential p, sT = p + p. It is apparent 
that the sudden decrease in p nearly equals the decrease in s, at fixed T and slightly changing constant p. If 
the temperature effects in the self- interaction potential V are neglected, ax ~ Jt — 0, then from Eq. we 
obtain the well-known relation between the critical temperature and the bag constant, B = (gq — g^) 



IV. DYNAMICS OF THE UNIVERSE DURING THE QUARK-HADRON PHASE TRANSITION 



The quantities to be traced through the quark-hadron phase transition are the energy density p, the tem- 
perature T and the scale factor a. These quantities are determined by the gravitational field Eqs. and ([7]) 
and by the equations of state ([15]), and ([T7| . We shall consider now the evolution of the Universe before, 
during and after the phase transition. 



A. Cosmological evolution in the quark phase (prior to the quark-hadron phase transition) 

Before the phase transition, at T > Tc, the Universe is likely in the partonic phase. With the use of the 
equations of state of the quark matter and of the Bianchi identity, Eq. ([7]) , the time evolution of the scale factor 
can be written in the form 

i^(T)^^ = -^^^|-i^|„ (23) 
a dttq 1 b aq 1 

and can be integrated to give the following scale factor-temperature relation: 



To' ' 
T 



- 1 



(24) 



where qq is the initial value of the scale factor corresponding to the temperature T = Tq of the universe, 
a (To) = flQ. In Fig. [U the variation of the scale factor of the Universe during the quark phase is presented as 
a function of the temperature T. Because of the expansion of the Universe the temperature is decreasing with 
the increase in the comoving time t. Therefore the scale factor a increases with the decreasing T. The exact 
numerical values of a{t) strongly depend on the initial Tq value. 

In order to have an analytical insight into the evolution of the cosmological quark matter, we consider the 
simple case in which the temperature corrections can be neglected in the self-interaction potential V. In this case 
V = B = constant and EoS of the quark matter is given by the bag model equation of state, Pq = [pq — 4i?) /3. 
Thus, Eq. ([7]) can immediately be integrated to give the following simple scale factor-temperature relation: 

Hence the presence of a temperature-dependent potential term V{T) in the quark matter EoS drastically modifies 
the scale factor-temperature relationship. The same result can be obtained by taking ax = 7t = in Eq. (j24p . 

With the use of Eq. p3)) and from the gravitational field equations, we obtain an expression describing the 
evolution of the temperature of the Universe in the quark phase, given by 

dT _ T3 ^(3a^-ar)T4 + 7rT2 + i3 ^^^^ 
dt V^mpi [(3ag - ar) /SaJ T^ + {'jr/Gaq) 

The variation of the temperature in the quark phase is presented, for different values of the bag constant B, in 
Fig. [21 The temperature dependence of the Hubble parameter during the quark-gluon phase is represented, for 
different values of the bag constant B, in Fig. [3| 
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140 160 180 200 220 240 

Fig. 1: The dependence of the scale factor a of the temperature T during the quark phase for Tq = 250 MeV. The a-T 
relation is almost independent on the mass nis of the strange quark. 
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Fig. 2: Time dependence of the temperature of the Universe T during the quark phase for a strange quark mass rUs = 200 
MeV and different values of the bag constant: B^/" = 100 MeV (solid curve), B^^^ = 200 MeV (dashed curve) and 
^1/4 _ MeV (long-dashed curve), respectively. 

B. Cosmological dynamics during the first-order quark- hadron phase transition 

During the quark-hadron phase transition, the temperature and the pressure are constants, T — Tc and 
p — pc, respectively. The entropy S = sa^ and the enthalpy W = {p + p) are conserved quantities. The 
energy density p (t) decreases from pq (Tc) = pq to ph (Tc) = pn- At the critical temperature Tc = 125 MeV, we 
have PQ ~ 5 X 10^ McV* and pH — 1-38 x 10^ MeV"', respectively. The value of the pressure of the cosmological 
fluid during the phase transition is pc ~ 4.6 x 10^ MeV^. Following [l^, it is convenient to replace p (t) by the 
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Fig. 3: Temperature dependence of the Hubble parameter H during the quark phase for a strange quark mass of 
ms = 200 MeV and different values of the bag constant: B^/" = 100 MeV (solid curve), B^'^ = 200 MeV (dashed curve) 
and B^''* — 250 MeV (long-dashed curve), respectively. 



volume fraction of matter in the hadron phase 

p{t) = pHh{t) +pq[1- h{t)] = PQ [1 + nh[t)] 



(27) 



where n = {pn — pq) I pq is the relative density and t is the comoving cosmological time. It is obvious that at the 
beginning of the quark-hadron phase transition, the quantity /i(tc) vanishes, where is the time corresponding 
to the beginning of the phase transition and p (tc) = pq- At the end of the quark-hadron transition, h (th) ~ 1, 
where is the time at which the phase transition ends corresponding to p{th) = pH- At t > th, the Universe 
enters in the hadronic phase. 

From Eq. ([7]), we obtain an expression for the Hubble parameter 



H 



1 eh 



1 {PH - Pq) h ^ 

^ pq + Pc + [pH - pq)h 3 1 + e/i' 



(28) 



where h denotes the time derivative of the hadron fraction parameter ft,, which can be utilized as an order 
parameter. 



^ ^ Ph - Pq 
Pq+Pc' 

Then, Eq. psp immediately leads to the scale factor, 

a{t) = a{Q{l+eh{t)y^/^ , 



(29) 



(30) 



where we have used the initial condition h (tc) = 0. The evolution of the fraction of the matter in the hadronic 
phase is described as 



1 



mpi 



3[l + nh{t)] pq 



h{t) + - 



with the general solution given by 

h{t) 



ne 



sech 



PQ 



t-t^ 
mpi 



(31) 



(32) 
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As given above, the quark-hadron phase transition ends up, when the value of h(t) reaches 1. Then, the time 
tfi at which the phase transition ends reads 

At the end of the phase transition the scale factor of the Universe has the value, Eq. ([5(1)) 

aith) = a{tc)ie + l)-^^\ (34) 

The variation of the hadron fraction given by Eq. p2p . as a function of the dimensionless time parameter 
X = y/7)Qtpit is represented, for different values of the parameter e and for n fixed in Fig. 21 The hadron 

fraction apparently gives an estimation for hadrons formed inside QGP. Having the expressions of h{t) and h{t), 
the analytical forms for both H and a can be directly obtained. It is straightforward to show that the Hubble 
parameter during the phase transition can be expressed as 

H = -jJ— J[l + nh{t)]pQ. (35) 
vim PI V 

The variation of the dimensionless Hubble parameter Hq — H mpi/^pg is represented, as a function of the 
dimensionless parameter x = ^/Pgtpity for different values of the parameter e and for n fixed in Fig. [5j 




X 



Fig. 4: Time evolution of the hadron fraction h during the quark-hadron phase transition for n — —0.74 and different 
values of e: e — —1/4 (solid curve), e — —1/2 (dashed curve), e = —3/4 (long-dashed curve) and e — —1 (very 
long-dashed curve), respectively. 



C. Cosmological evolution in the hadronic phase (post quark-hadron phase transition) 

Finally, after the phase transition, the energy density of the pure hadronic matter is ph = iph = SattT'*. The 
Bianchi identity Eq. (O gives 

a{T) = a{t,,) ^. (36) 



11 




Fig. 5; Time evolution of the dimensionless Hubble parameter Ho = H x mpi/^/pQ during the quark-hadron phase 
transition as a function of the dimensionless time parameter x = ^PQ^Pit for n = —0.74 and different values of e: 
e — —1/4 (solid curve), e — —1/2 (dashed curve), e = —3/4 (long-dashed curve) and e = —1 (very long-dashed curve), 
respectively. 



The time evolution of the temperature in the hadronic phase is governed by the equation 

§ = (3a. T^f' = - — V^T\ (37) 

at y/3mpi mpi 

giving a comoving time 

i — i,- = 

2 ^ V T2 
From Eq. and (P7)) . the Hubble parameter reads 



t-t.-^{^^-l\. (38) 



H[T) = Y:^t\ (39) 
mpi 

H{t) = ^ i ^ . (40) 

mpi + {2^lmpi) {t - th) 



During the hadronic phase, the density of the Universe varies wit\\ the time as 



P>^ii) = ^—^, -72- (41) 



[T^^ + {2^/mpi) {t-th)\ 



The temperature dependence of the scale factor a of the Universe during the hadronic evolution phase is 
presented in Fig. [6l The temperature dependence of the Hubble parameter H is represented in Fig. [71 

Finally, in Fig. |5]we present the time evolution of the scale factor a of the Universe during the quark phase, the 
phase transition and the hadron phase, respectively, for several values of the bag constant B. We assume that 
the quark phase begins at a time t = tq, when the value of the scale factor of the Universe is a = a {tq). The 
phase transition temperature is assumed to be Tc = 125 MeV, with a corresponding quark matter energy density 
at the transition moment of /5q = 5 x 10^ MeV''. For the parameter e we have taken a value of e = — 1/4. As 
one can see from the figure, an increasing value of the bag constant accelerates, in the long term, the expansion 
of the universe. 
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Fig. 6: The T dependence of the scale factor a during the hadronic phase. 
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Fig. 7: The T-dependence of the Hubble parameter H during the hadronic phase. 



V. PHASE TRANSITION IN LATTICE QCD SIMULATIONS AND HEAVY-ION COLLISIONS 

Before introducing the QCD EoS, it is useful to study the similarities between heavy-ion collisions and the 
early Universe Q . It is conjectured that the first-order phase transition, studied in scction lllll might take place in 
the heavy-ion collisions and/or in lattice QCD simulations. Such a prompt transition seems to have fundamental 
astrophysical consequences. Its dynamics has been discussed in the previous section. Despite of the order of 
phase transitions, the QGP era seems not to be followed by an extreme expansion (inflation). This is apparently 
the case in heavy-ion collisions, because of the baryon number conservation and the limitation of baryon-to- 



photon ratio (n;, — ni)/n^ 



IQ- 



[481. Therefore, np 



recently measured by ALICE experiment at 7 GeV 



can be used to estimate the photon number density, ~ 5.5 x 10 , while in the CMB era, ~ 411.4(T/2.73K) 
cm~^. Furthermore, the QGP era seems to be the last symmetry-breaking era of strongly interacting matter. 
By symmetry breaking, we mean deconfinemcnt and chiral symmetry breaking and/or restoring, respectively. 
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Fig. 8: Time evolution of the scale factor during the quark phase (solid line, quark-hadron phase transition (dashed line) 
and hadron phase (long-dashed line), respectively, for several values of the bag constant B: B^^^ = 100 MeV, B^^"^ = 200 
MeV and B^^* = 250 MeV. The numerical values of the scale factors raise with increasing B. The quark phase begins 
at t = tq = 0, when the value of the scale factor is given as ao = a (tg) = a(0) — 10^*. The assumed critical transition 
temperature is Tc — 125 MeV, the quark density is pg = 5 x 10^ MeV*, while e — —l/A. 




In an isotropic and homogeneous background, the volume of the Universe is directly related to the scale factor 
a{t), where t is the comoving time. Implementing a barotropic EoS for the background matter makes it possible 
to calculate - among others - the Hubble parameter H{t) = a{t)/a{t). Focusing the discussion on QCD era of 
the early universe, which likely turns to be fairly accessible in high-energy experiments, the equation of state is 
very well defined. In Ref. [40], a viscous EoS for QGP matter has been introduced and different solutions for 
the evolution equation of H have been worked out. The ratio of baryon density asymmetry to photon density, 
ry, has been measured in WMAP data Then Up — Up from the ALICE experiment at 7GeV can be used 
in order to estimate the photon number density, ~ 5.5 x IC*, while in the CMB era, ~ 411.4(r/2.73K) 
cm~^. 

In a comoving volume V ~ a'^(t), the number density of noninteracting photons is supposed to remain 
constant. Therefore, ~ l/a^{t). Nevertheless, when the Universe was expanding, T decreases and a^{t) n-y has 
to be affected. The previous values of support this conclusion. There is a conserved quantity accompanying 
such a transition, namely, the entropy density s. In a perfectly closed system like the universe, s likely remains 
unchanged. From the first-law of thermodynamics [40| one can show that at vanishing chemical potential, 

, , P(T) + p(T) ai + 1 , „ 
1 a2 

implying that s{T) is related to p(T)l-'^^ where oi = 0.319, 02 = 0.718 ± 0.054 and 03 = 0.23 ± 0.196. This 
relation is valid at low energy, where the dominant degrees of freedom are given by hadron resonances. Baryon 
and boson relative abundances {n{T) — n{T))/{n{T) + n(T)) can be studied in hadron resonance gas (HRG) 
model. It is found that the abundance approaches 10~^. For instance, the kaon relative abundance is by about 
one order of magnitude higher than that of the proton. It is obvious that the abundances of the li ght elements 
(^Li, ^He, '^He and ^H) produced in the early Universe are sensitive indicators of number density [43|. Recent 
lattice QCD calculations [13, [EH give estimations for the equation of state, temperature and bulk viscosity of 
hadronic and partonic matter at high temperatures. 

As we will show below, the gravitational cosmological field equations, Eqs. (|4]) and ([7]), relate the cosmological 
parameters, like the Hubble parameter H and the scale factor a, to the energy density p. Again, the barotropic 
equations of state for the thermodynamic parameters are standard in analyzing the viscous cosmological models, 
whereas the equation for r is a simple procedure to ensure that the speed of viscous pulses does not exceed the 
speed of light. 
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Using this set of equations seems to define the vahdity of this treatment. Apparently, it depends on the 
validity of the equations of states, which we have derived from the lattice QCD simulations at temperatures 
larger than Tc — 0.19 GeV. Below Tc, as the Universe is cooled down, not only the degrees of freedom suddenly 
increase [s^ . but also the equations of state turn out to be the ones characterizing the hadronic matter. Such a 
phase transition - from QGP to hadronic matter - would characterize one end of the validity of our treatment. 
The other limitation is provided by the very high temperatures (energies), at which the strong coupling as 
entirely vanishes. 

The lattice QCD simulations benefit from the rapid progresses achieved in computational facilities and al- 
gorithms. The accuracy of recent lattice results is comparable with the laboratory experiments. Currently, 
it is possible to perform lattice QCD simulations at almost physical quark masses. Recent results on QCD 
equations of states have been reported in (ssj . It is apparent, that the influence of radiation and leptons on 
phase transition are minimized 54 1 . 



A. Equations of state of the viscous quark-gluon plasma 

In this section, we give a list of barotropic equations of states deduced from the lattice QCD simulations [H^l 
(an analytic crossover phase transition is obtained) and the quasiparticlc model [s^ . The latter is utilized when 
no lattice QCD results are available. Figure [S] depicts the pressure p in dependence on the energy density p in 
a wide range of temperatures, 1/2 < T/Tc < 3. Details on the lattice configurations are described in [s^. It is 
obvious that the (barotropic) pressure - energy density dependence is almost linear referring to the nature of 
the phase transition from hadrons to quarks and vice versa. This confinement-deconfimement phase transition 
seems to be smooth i.e., simply continuous and takes place very slowly. This kind of transitions is a very 
moderate than the second-order one. The nature of the phase diagram in lattice QCD has been discussed 
in (56| . In Fig. ^ the dashed line represents the fitting in the entire T-region. Ignoring the dip around Tc, the 
results can be fitted as a power law, 

p(p) = alp"^ (43) 

where ai = 0.178 ± 0.009 and a2 = 1.119 ± 0.011, respectively. In the hadronic phase i.e., at temperatures 
< Tc, the previous power law dependence seems to remain valid. Little changes appear in the parameters; 
ai = 0.096 ± 0.003 and a2 ~ 1.03 ± 0.04. In the quark phase i.e., at temperatures > Tc, the following 
polynomial 

p(p) = -^ + aip"% (44) 

describes this barotropic equation of state, where ai = 0.221 ± 0.004 and a2 = 1.072 ± 0.005. 

Apparently, expressions ([^5)) and imply that the speed of sound drastically changes with the changes in 
the phases: 

• hadron and quark phase: = dp /dp ~ 0.199 pO^^ i.e., the sin the peed of sound depends on the energy 
density. The latter has a nonmonotonic behavior when going from hadronic to partonic phases and vice 
versa. 

• in the hadron phase: ~ 0.098, 

• in the quark phase: ~ 0.237. 

Figure [TUl presents the barotropic dependence of T as calculated in lattice QCD. The relation can nicely be 
fitted by the polynomial 

T(p) = /3i + /?2/^ (45) 

where /3i = 0.123 T 0.004, ^2 = 0.058 ± 0.0038 and /Jg = 0.39 ± 0.013. Again, the specific heat in the whole 
phase of hadrons and quarks seems to depend on the energy density, strongly 
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Fig. 9: The pressure density, p, is drawn in dependence on the energy density p. Both quantities are given in physical 
units. Symbols are lattice QCD calculations using p4 action and temporal lattice size, Nt = 8 [S^l- The dotted curve 
is the fitting in the quark phase, Eq. (|44|l . The dash-dotted curve gives the fitting in the hadronic phase. The overall 
fitting is given by Eq. (|43|l (double-dotted curve). The small dip at Tc seems to reflect the slow phase transition known 

as crossover . 
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Fig. 10: The symbols are the lattice QCD calculations for the barotropic dependence of temperature T{p) using p4 
action and temporal lattice size, A''r — 8. Details on the lattice configurations are given in Ref. [s^l. The solid curve 
gives the fitting according to Eq. (|45)l . 



In determining this value, the volume V is conjectured to remain unchanged. 

Based on the quasiparticle model [l^l , which is an effective model used to reproduce the lattice QCD results 
on various thermodynamic and transport properties (ssj . the bulk viscosity ^ reads 



ir_ _ f _ rr,de \ dp 
3e V dTj dp 



(47) 
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where g is the degeneracy factor of quarks, giuons and their antiparticles. The function /o gives the distribution 
function (exp(e) ± 1)~^ of boson and fermion particles, respectively. The quantity t—[p^-\- n(T)]^/^ stands for 
the effective dispersion relation of single particles. It depends on the particle mass n(r) which in turn varies 
with the effective coupling G{T). Therefore, the effective coupling G{T) plays an essential role in this model. 
It has to be adjusted to reproduce lattice QCD results. In top panel of Fig. [Til the bulk viscosity coefficient ^ 
is given as a function of T. The results are fitted very well as 

= In (^^) [73 + 74 (T - 75)"1 . (48) 

At high T, the fitting parameters are 71 = 1.042 ± 0.067, 72 = -0.0035 ± 0.0059, 73 -0.329 ± 0.058, 
74 = 25.666 ± 1.521 and 75 = -0.367 ± 0.0159. At low T, the fitting parameters read 71 = 0.801 ± 0.595, 
72 = -0.352 ± 0.558, 73 = -0.350 ± 0.174, 74 = 303.582 ± 3.126 and 75 = 0.189 ± 0.054. In this region the 
exponent 76 = 1. In the high-T region, 7g = 6. 

Again, in quasiparticle model j55l l57l|. the relaxation time reads 

.-.,r,.^TO'(r),„(j^). ,49) 

where 0,^ = 6.8 [s^l- The results arc drawn in bottom panel of Fig. [TTJ Fitting of these results leads to 

where Si = 2.362 ± 1.318, ^2 = -0.022 ± 0.056, ^3 = -3.176 ± 1.05, 64, = 0.435 ± 0.126, 65 = 2.362 ± 0.318, 
(56 = 3 and 67 = 1.25 ±0.12. 

From the two expressions (|48|) and (f50|) . it is obvious that the barotropic relations of ^ and r are related to 
each other. Such a relation has been modeled by the projection operator method [H^ as 

Furthermore, the bulk stress is to be related to the distribution function of relaxation time (60j . Such a 
dependence has to follow the causality principle and fits perfectly with the laws of thermodynamics [6l|. The 
speed of sound = dp/ dp, can be taken from the lattice QCD simulations [5^. The results of c^(T) are given 
in Fig. [121 Below Tc, the lattice results show a small peak. Remarkable work has been devoted to accurate 
its location and altitude. The results from the hadron resonance gas model are given, as well. Although the 
appearance of the peak, the disagreement is not to be neglected. With reference to the restricted causality 
principle, the nonmonotonic behavior of below Tc would be explained in the light of: 

• baryon and strange degrees of freedom which would play an essential role in reproducing c^(T, ^j,), where 
fif, is the baryo-chcmical potential, 

• the interpolation of both entropy .s(r, fit,) and specific heat cv{T, fib) which has been suggested to partly 
explain nonmonotonic behavior below Tc, 

• the condition(s) deriving the chemical and thermal freeze-out which would enlighten such a behaviour, 

• the interactions between the constituents of the hadronic phase are conjectured which would be able to 
explain the nonmonotonic entropy and specific heat production and 

• the time- varying equation of state in the hadronic phase which refers to out-of-equilibrium processes, while 
their modification in thermal and dense matter would refer to symmetry changes. 



B. Bulk Viscosity in the Hadronic Phase 



The treatment of bulk viscosity in Hagedorn fluid has been studied in Rcf. [1^. Such a fluid is conjectured 
to be composed of hadrons and resonances with masses m < 2 GeV. The treatment is based on the relativistic 
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Fig. 11: The top panel depicts the dependence of the bulk viscosity on the temperature T as deduced from the 

quasiparticle model (symbols) [H, [sj- The dotted curve gives the fitting of results related to the hadronic phase. The 
fitting of the results in the partonic phase is given by the dashed curve, Eq. (|48|l . The two regions meet at Tc ~ 270 MeV. 
Amazingly, this value finds its root in lattice QCD. The quenched lattice QCD calculations, where the quark masses are 
supposed to be very heavy, predict that the critical temperature separating hadrons from QGP has the same value. In 
bottom panel, relaxation times r is drawn against T in MeV units. The solid curve gives the fitting according to Eq. 
(|50|l . Again, there are two separate regions. In the first one, r raises with increasing T. While, in the second region, r 
decreases with increasing T. Such a nonmonotonic behavior is characterized at Tc ~ 270 MeV. 

kinetic theory formulated under the relaxation time approximation. The in-medium thermal effects on bulk (and 
shear) viscosity and the averaged relaxation time with and without the excludcd-volumc approach are deduced. 
It has been suggested that the dynamics of the heavy-ion collisions, the noncquilibrium thermodynamics and 
the cosmological models, require thermo- and hydrodynamic cquation(s) of state. When assuming vanishing 
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Fig. 12: The speed of sound = dp/dp is drawn versus T in the physical units. The symbols represent the results from 
lattice QCD calculations [H^l, using p4 action and temporal size A'^t ~ 8. Below Tc, results from the hadron resonance 
gas model are given (solid curve on the left). The other two curves are fittings for lattice QCD data. 



chemical potential and the heat conductivity, the bulk viscosity in thermal medium reads 

m = ^^^p(m,)^°°no(l + no)(^c^e,?-ip2)^^dp, (52) 

where p{mi) is the Hagcdorn mass spectrum p{m), which implies growth of the hadron mass spectrum with 
increasing the resonance masses. 

p(m) = A (toq + TO^)'^^^ exp(m/T/f ), (53) 

with k ^ -5, A ^ 0.5 GeV^/^, mo = 0.5 GeV and Th = 0.195 GeV. The number density no is related to the 
deviation of energy-momentum tensor from its local equilibrium ST'^'^ . Such a deviation is corresponding to the 
difference between the distribution function near and at equilibrium, 6n = n — no- The latter can be determined 
by the relaxation time approximation with vanishing external and self-consistent forces [62l . l63j 

6nip,T) = -T{T)^d^no{p,T). (54) 
p ■ u 

The nonequilibrium number density n{p,T) is to be decomposed using the relaxation time approach into n = 
no -I- TUi -I- • • • . Alternatively, as n{p, T) embeds the l**-rank tensor u, 5T^" can be decomposed into u [g^ in 
order to deduce its spatial components. 

The relaxation time depends on the relative cross section as 

"(^^ = nf{T){v{T)a{T)y ^''^ 

where v(T) and nf(T) is the relative velocity of two particles in case of binary collision and the density of each 
of the two species, respectively. The thermal-averaged transport rate or cross section is {v{T)a{T)). 

The ratio of bulk to shear viscosity, ^/tj, can be related to the speed of sound in a gas composed of massless 
pions. Apparently, there are essential differences between this system and the Hagedorn fluid. According to 
[64| . the ratio of £,/r] in N ^ 2* plasma is conjectured to remain finite across the second-order phase transition. 
In Hagedorn fluid, the system is assumed to be drifted away from equilibrium and it should relax after a 
characteristic time r. Should we implement a phase transition in Hagedorn fluid, then r cx where z is the 
critical exponents, which likely diverges near Tc- 
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C. Deconfinement and chiral phase transitions (crossover ) in lattice QCD simulations 

Remarkable advances have been made in studying the equilibrium properties of the phase transitions. Obvi- 
ously, the phase transition is coupled with symmetry breaking and out-of-equilibrium. Therefore, it is natural 
to turn our attention to the consequences when the system is enforced to go through an out-of-equilibrium phase 
transition. Thermodynamically, the first and second order phase transitions are described by continuous first 
and second derivative of the free energy, respectively. The infinite order phase transition is also continuous. But 
it breaks no symmetry. A famous example for it is the Kosterlitz-Thoulcss transition in the two-dimensional 
XY model [63 ■ The crossover phase transition of lattice QCD simulations is likely a continuous one. An 
out-of-equilibrium state is reached when the system in deviated from its equilibrium state by applying an in- 
stantaneous perturbation. The system will relax to its equilibrium state by dissipating the energy transferred 
during the transition j66| . Relating the amplitude of this dissipation to the amplitude of fluctuations in equilib- 
rium dates back to Einstein's work on the Brownian motion in 1905. Lars Onsager established a hypothesis that 
the relaxation of a macroscopic nonequilibrium perturbation follows the same laws which govern the dynamics 
of fluctuations in equilibrium systems. In other words, the regression of microscopic thermal fluctuations at 
equilibrium follows the macroscopic law of relaxation of small nonequilibrium disturbances [67j . 

The lattice QCD simulations turn out to be an accurate tool to study - among others - the thermodynamics 
of the hadronic and partonic matter up to temperatures of couple Tc] the critical temperature Tc [26l - l3^ . For 
two quark flavors (n/ = 2) the phase transition is second-order or a rapid crossover . Tc ~ 173 ± 8 MeV. 
For n/ = 3, the phase transition is first-order and Tc ~ 154 ± 8 MeV. For nj = 2 + 1, the transition is again 
crossover and Tc ~ 173 ± 8 MeV. For the pure gauge theory, Tc ~ 271 ± 2 MeV and the phase transition is 
first-order. In all these lattice QCD simulations, the quark masses are heavier than their physical values. At 
physical masses, the critical temperature for ri/ = 2 -I- 1 is ~ 200 MeV. Apparently, we conclude that the order 
of the deconfinement phase transition can be either first, or second or crossover (infinite). It depends on the 
quark flavors and their masses. The extreme conditions in the early Universe likely affect the properties of the 
hadronic and partonic matter. 

The chiral phase transition is assumed to accompany the deconfinement one, especially at vanishing chemical 
potential. It is expected that the restoration of the chiral symmetry breaking takes place in full-perturbative 
and nonperturbative QCD at high temperatures, if the matter is assumed to be exclusively built of light and 
strange quarks (6^. In perturbative QCD, the chiral symmetry is valid for massless quarks. It is entirely broken 
in the hadronic phase. It not yet completely clear what is the order of phase transition between hadronic and 
partonic QCD phases when the broken symmetry is restored at finite temperatures and densities. Different 
lattice QCD simulations, mainly referring to chiral condensate and chiral susceptibilities j69| indicate that the 
chiral phase transition is of the second order at vanishing chemical potential: 

S\J{nf)i X S\J{nf)r ^ S\J{nf )v- (56) 

The chiral condensate vanishes at the limit mq —J- [t^. Below Tc, the chiral condensate entirely vanishes, as 
well. It is finite above Tc, 

{i^^)=-^J-lnZ, (57) 
V oruq 

where In Z is the partition function describing the system. The chiral perturbation theory proved to be a very 
important method in determining some essential observables in QCD, which are dominated at low temperature, 
such as the masses of pseudoscalar mesons, their decay constants and the chiral observables. It provides an 
explanation why pseudoscalar mesons are very light. The Goldstone theorem states that for each generator of 
a spontaneously broken symmetry, there exists a massless Goldstone boson with spin and with symmetry 
properties that are related to those of the symmetry transformation. The Goldstone bosons of the chiral 
perturbation theory are just the pseudoscalar mesons. This can be utilized as a signature for the phase transition. 

VI. DYNAMICS OF THE BULK VISCOUS QUARK-GLUON PLASMA FILLED UNIVERSE 

In the following we consider the cosmological evolution of the viscous quark-gluon plasma filling the Universe 
in the framework of both Eckart and the full causal approaches to dissipativc processes. 



20 



A. Evolution of the Hubble parameter in the Eckart model 



Substituting the barotropic expressions (|43p . (|45p and (|48p in Eq. ([TT|) and by assuming for the bulk viscous 
pressure the Eckart relation, given by Eq. (|10p , we obtain for the evolution of the Hubble parameter the equation 



2 A I 



3 In 



132 A^^ ij2/93 



71 



72 



73 + 74 (/32 H^^' - 75) 



76 



0, (58) 



where A = 3/(87rG'). As given in the introduction, the Planck units are given by this parameter, A = 1.778x 10' 
GeV^. This differential equation can be solved, analytically, when assuming that 



,37 



i/'^^ =(7l+72)//32A^^ 



(59) 



Then, in terms of H, the comoving time reads 



t = 2A 



In [-A (1 + aa - 3/3274) + 3 (73 - 7475) /H] 
9(73 - 7475) 



(60) 



Figure [T3l shows the dependence of t on H. It describes a universe, where its background fluid is characterized 
by Eckart theory. The three curves differentiate between different types of the background matter. A discussion 
about the effect of background matter is given in section IVII Al The collisionfree and nonviscous background 
matter is given by the dashed curve. Solid and dotted curves describe the t — H relation when the background 
geometry is filled with viscous hadron-QGP and nonviscous QGP, respectively. At small H values, there are 
obvious differences between the latter types of matter and between them and the ideal matter. At large H 
values, the comoving time behaves very smooth with H , although hadron-QGP results in larger t than QGP. 
In both of them, t is larger than in ideal matter. 
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Fig. 13: The comoving time t is drawn (solid curve) in dependence on the Hubble parameter H, Eq. (|60p . The 
background fluid is characterized by Eckart theory. The dotted curve gives the results when viscous QGP equation of 
state is implemented [3343^ . ItH . It^ ] . The dashed curve shows the results when the background geometry is assumed to 
be filled with an ideal gas. The Planck scale is given in physics units. 
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B. Evolution of the Hubble parameter in the full causal approach 



Comparing Eq. ([5T|) with the expressions P5|) and ([50)) makes it quite apparent to have a barotropic expression 
for the relaxation time [4l| 




Fig. 14: Parametric functions given in Eqs. H63p and H64p are depicted in dependence on H. This illustration gives an 
indication about the dependence of the function g on the independent functional parameter z. Such a dependence is 
given in Fig. (fTS]) . 
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Fig. 15: The parametric dependence of the function g on the independent functional parameter z. Apparently, it has a 
nonmonotonic behaviour. At _ff-values, where this study is valid, the dependence is apparently linear, g{z) oc z. 



Witli the use of Eqs. (|43p . pS)) and p^ . respectively, we obtain the following equation describing the 

cosmological evolution of the Hubble parameter H: 




7 1 3ai 

2 ~ 27r ^ ItT 



_3_ 

2^ 



473 + /32j4\ - H - 47475 In 



47r if 4^ 47r 



473 + I32l4\ - H ~ 47475 In 




y HH- 



= 0. (62) 



In obtaining this equation we have introduced a number of very tiny approximations. To the exponents a2^ ^3 
and 76 we have assigned the values, 1, 1/2 and 1, respectively. In order to derive an analytical solution for this 
Abel differential equation, we follow the procedure given in |37l . [7l| . After some Algebra, it ends up with these 
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two functions 

9iH) = 



14 7r-2-3ai+47r/| 473 + /?2 74 V^A-H^ - 47475 In {ji2^/&fn H - A-flj / 'i-f2 |' 



z{H) = 



which are plotted in Fig. p4)) . They play an essential role in deriving an analytical solution for Eq. (|T4| 
Approximating the parametric dependence of g{H) on z{H), Fig. [TSl wc get the linear dependence, 



(47r-l)(87r-l)/3|7| 

|47r(l - 87r)73 + 47r(87r - 1) [7475 + (73 - 7475)] + /3274(47r - 1) Vg^ + {Utt - 2 - 3ai)^274 



Then, from the definition of fl, we simply derive 



0.192 z. 



(65) 



(66) 



In order to reduce this expression to the canonical equation of Abel type, we use the relation = z/V. Then 
from Eqs. and ([M)) . we obtain a first-order differential equation for H, 



VH = 



(47r-l)(8^-l)/3|7l 



H {47r(l - 87r)73 + 47r(87r - 1) [7475 + (73 - 7475)] + /3274(47r - 1) + (147r - 2 - 3ai)/3274] H] 

where P is a free parameter. The solution simply reads 

(87r - l)/?274 V 



t = 



/6^ + (14 TT - 3ai - 2) /3274] H 



(68) 



The dependence of the cosmological comoving time t on the Hubble parameter H is graphically illustrated in 
Fig. 1161 It is apparent that t{H) is monotonic. The same dependence has been obtained, when assuming that 
the background matter is characterized as an ideal gas, t = 2/{3'-fH). All this is summarized in Fig. [161 Solid, 
dashed and dotted curves represent the results for viscous hadron-QGP, viscous QGP and ideal (nonviscous 
and collisionless) matter, respectively. 



VII. COSMOLOGICAL IMPLICATIONS 



Assuming that the background geometry is filled with Eckart relativistic viscous fluid, the comoving time t is 
given as a function of the Hubble parameter H in Eq. (|60| and drawn in Fig. [13] (solid curve). The dotted curve 
gives the results when viscous QGP EoS is implemented j33 - l38l . ItiI It^. The dashed curve draws the results 
when the background geometry is assumed to be filled with an ideal gas. In Fig. [1^1 another t-H dependence 
is obtained when assuming that the cosmological background is filled with Israel-Stewart relativistic viscous 
fluid. The solid curve represents the results of present work, in which the background matter is assumed to be 
characterized by viscous hadrons and QGP i.e., including phase transition(s). The same treatment is applied 
for viscous QGP and ideal gas. The results arc drawn by dotted and dashed curves, respectively. 

Before discussing the potential cosmological implications, it is in order now to elaborate essential aspects. 
We start with the phase transition in the early universe. The flrst-order phase transition has been discussed in 
section IIIII Section [Vj was devoted to discuss the phase transition(s) as measured in lattice QCD simulations. 
Accordingly, we conclude that the order of the conflnement-deconflnement phase transition depends among 
others on the effective degrees of freedom and the matter content (quark flavors, etc.) The results given in 
Figs. [T3| and [T6| illustrate the effects of degrees of freedom (ideal gas, QGP and hadron-QGP matter) and 
in indirect way the phase transition (QGP matter above Tc and hadron-QGP matter over a wide range of 
temperatures). The evolution of the Hubble parameter obviously depends on all these factors. This might 
have a direct cosmological implication that our picture about the expansion of the Universe has to be revised, 
accordingly. 
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Fig. 16: The cosmological comoving time t vs Hubble parameter H, Eq. (|68|l . is graphically illustrated. The treatment 
of background matter is done by Israel-Stewart theory. The solid curve represents the results of present work, where the 
background matter is assumed to be characterized by viscous hadron and QGP. The dotted and dashed curves give the 
results when viscous QGP and an ideal (coUisionless and nonviscous) gas, respectively, are assumed to fill the background 
geometry. Planck scale is given in physical units. 



Other cosmological implications might arise as a consequence of the phase transition itself. The first-order 
phase transition is to be characterized by a sudden change in the symmetry. It exhibits a discontinuity in the 
first derivative of the free energy with respect to some thermodynamic variable. In the cosmological context, 
such a transition is accompanied by bubble nucleation [tI]. In light of this, the Universe is conjectured to go 
from a metastable state to a new phase, a true vacuum state through the nucleation of bubbles of the new 
state (rsj . Implementing this model to the hadron-QGP transition makes it possible to suggest a scenario, in 
which the Universe starts from QGP state and ends up in the hadronic state through the nucleation of hadrons. 
Depending on the kinematics of the bubble nucleation, the Universe might or might not "recover" from this 
type of phase transition and its relics are left behind i.e., relic QGP objects. The latter would survive for a very 
long time. The abundance and the size of the quark nuggets have been discussed in [2l| . Objects with a quark 
content ranging from 10~^ to 10 A/q could have been formed during the cosmological phase transition. 

Furthermore, a significant amount of entropy production is to be released during such a process, so that 
at vanishing chemical potential s = (c^ -f 1) p/T. The density fluctuations are assumed to be amplified by 
vanishing speed of sound during the quark-hadron phase transition. Fig. 1121 The lattice QCD and hadron 
resonance gas calculations show that the speed of sound reaches a minimum value, ~ 0.1, at Tc- On the 
other hand, the density fluctuations could produce QGP lumps decoupled from the expansion, which rapidly 
transform into quark nuggets. Typical distance between bubble centers is conjectured to be of the order of a few 
meters. It is worthwhile to mention here that the resulting baryon inhomogeneities may affect the primordial 
nucleosynthesis. Such a cosmological consequence can be observed. The origin of inhomogeneities in the matter 
distribution, which are assumed to be responsible for the later formation of galaxies, cannot be explained by 
density fluctuations, alone. After flxing the baryon number, the appearance of these fluctuations is almost purely 
adiabatic. Any departure from adiabaticity falling off is inversely proportional to the mass of the perturbation 
(t^ . This will be elaborated in next paragraph. 

At the phase transition, the scale of the cosmological QCD transition is assumed to be given by the Hubble 
radius Rh- Quantitatively, Rh — mpi/T^ ~ 10 km. The mass inside the Hubble volume is ~ IM©. At the 
QCD phase transition, the expansion time scale is 10~^ s, which is much large in comparison with the time-scale 
of QCD, 1 fm/c ~ 10~^^s. Even the rate of weak interactions seems to exceed the Hubble rate by a factor 
of lO''. Therefore, we conclude that photons (radiation), leptons, quarks (fermions) and gluons (bosons) are 
lightly coupled and may be described by an adiabatically expanding fluid [40l [72j , as the transition takes place 
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in an extremely short time. 

The current heavy-ion experiments program, LHC, seems to be very close to probe early eras of the universe. 
It seems to produce similar antiparticle and particle, when not entirely identical [TSj . This can be taken as 
another supportive indicator for utilizing EoS deduced from heavy-ion collisions and/or lattice QCD calculations. 
It seems that the observed matter-antimatter asymmetry can be explained without recourse to the hypothesis 
of specific initial conditions fl^ . 



A. Different types of background matter 



The different phase transition(s) likely change the symmetries and thereupon different phases or types of 
matter are to be expected. The dynamics of the Universe during the fist-order phase transition from QGP to 
hadrons has been discussed in Section lTVl Such a transition is assumed to go over three phases defined by various 
symmetries. Prior to the phase transition i.e., partonic (QGP) symmetry, the evolution of some cosmological 
parameters (H, a and p) have been studied by using the Bianchi identity. To have an analytical insight into 
the evolution, T-corrections are neglected in the self-interaction potential. The second phase deals with the 
dynamics of the Universe during the phase transition i.e., mixed phase symmetry. Here, T and p are assumed 
to remain unchanged. The entropy s and enthalpy W remain conserved, as well. The third phase is the one 
in which the dynamics of the Universe is studied post quark-hadron phase transition era (hadronic symmetry) . 
First we start with the time evolution of T. Then, we estimate the comoving time t. Again, the Bianchi identity 
helps us in expressing the scale factor a and the Hubble parameter H. The time evolution of the hadron fraction 
h describes the conversion process of QGP into hadrons. Therefore, it can be taken as a parameter describing 
the phase transition itself. 

Again, in this type of transition (a first-order phase transition through hadron nucleation), the numerical 
estimation of the cosmological parameters gives a clear indication that their time evolution varies from phase to 
another. In the QGP phase, the scale factor a normalized to ag is much smaller than that in the hadronic phase, 
(compare Fig. [T] with the top panel in Fig. IH]). When studying the Hubble parameter H, the T-dependence is 
just the opposite of a(T), (compare Fig. |3]with the bottom panel in Fig. 

This behavior can be compared with the case of another types of phase transitions, crossover . In Figs. [T3] 
and [ini we notice that the time evolution of H also depends on the type of matter filling the background 
geometry. If it is filled with QGP, the values of H are relatively large. It is relatively small, if the background 
geometry is filled with quarks and hadrons, especially when crossover phase transition is allowed to take place. 
Consequently, it is likely to predict that H in the hadron era is smaller than its value in the previous eras: 
mixed phase of partons hadrons and QGP. 

It seems to be in order now to highlight the differences between viscous and nonviscous background matter. 
By eliminating the dynamics controlling the phase transition, for instance, we assume that the background 
geometry is only filled with QGP. For simplicity, we utilize the Eckart theory. A comparison is illustrated in 
Fig. [T71 We notice that the viscosity seems to drastically slow down the evolution of the Hubble parameter. 
Should this result be confirmed, it would mean that the whole picture about the evolution of early Universe has 
to be revised. As a prompt consequence, one would expect a considerable delay in all phases post to QGP era. 
In order to make an estimation for this effect, other initial conditions have to be taken into consideration, for 
example, dynamics of phase transition(s), interaction(s), out-of-cquilibrium processes, etc. 



VIII. DISCUSSIONS AND FINAL REMARKS 



In natural units, ft = c = fc^ = 1, all expressions are given the Planck mass nipi. We consider the cosmic 
evolution of the early Universe in the regime of confinement QCD phase transition taking finite bulk viscous 
effects into account. Thereby, it is assumed that the bulk thermodynamic quantities are dominated by the 
strongly interacting matter component. Two cases, a first-order phase transition scenario and an analytic 
crossover transition, are considered. In this respect, the present work continues a previous series (33 - l40l . [tiI [72j 
in several aspects. Refined equation(s) of state based on newer lattice QCD results are considered. Different 
bulk viscosity expressions based on quasiparticle model are used. Finite cosmological constant has been utilized 
in Ref. [s^. Moreover, the infiuence that a first-order phase transition (neglecting viscous effects) is elaborated 
in the present work. 

Many details of QCD phase transition(s) are not yet conclusively understood. Even the order of transition 
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Fig. 17: The t — H relation in Eckart relativistic QGP fluid. Solid and dotted curves represent viscous and nonviscous 
QGP fluid, respectively. It is obvious that the viscosity slows down the evolution of the Hubble parameter. 



is still a matter of debate. An advance in understanding the numerieal values of the QCD coupling constants 
would be very helpful in obtaining accurate cosmological conclusions 38[. Such an advance may also provide a 
powerful method for testing on a cosmological scale the theoretical predictions of the brane world models and 
the possible existence of the extra-dimensions. Furthermore, the critical temperature Tc has been a subject 
of different lattice QCD calculations [26l - [3^ . In addition to this, it is still an open question whether both 
deconfinement and chiral phase transitions take place at the same Tc. 

The cosmological behavior in first-order phase transition can be characterized as follows. At the critical 
temperature, the energy p and entropy s densities decrease, suddenly. At fixed T and constant p, both quantities 
have the same rate. Depending on the symmetries, the transition is assumed to go through three phases. Prior 
to the phase transition i.e., partonic (QGP) symmetry, the evolution of some cosmological parameters {H , a and 
p) have been studied by using Bianchi identity. To have an analytical insight into the evolution, T corrections 
are neglected in the self- interaction potential. The second phase is the one during the phase transition i.e., 
mixed phase symmetry. Here, T and p are assumed to remain unchanged. Also the entropy s and enthalpy W 
remain conserved. The third phase is the one in which the dynamics of the Universe is studied post quark-hadron 
phase transition era i.e., hadronic symmetry. The Bianchi identity helps in expressing scale factor a and Hubble 
parameter H . The behavior of a and H with the cosmological comoving times follows the standard cosmological 
model. Both quantities are expressed in terms of the fraction of matter. The latter gives an estimation for 
hadrons that are formed inside QGP. The time evolution of the hadron fraction describes the conversion process 
of QGP into hadrons. Therefore, it can be taken as a parameter describing the phase transition itself. A 
quantitative comparison between the evolution of scale factor a in the three phases show that a increases while 
moving from quarks to hadrons over the mixed phase. The values of the bag pressure are reflected in these 
calculations. In all phases we find that increasing the bag pressure raises the value of the scale factor. 

Taking into account the recent lattice QCD results, we find that the order of the phase transition can be cither 
continuous or discontinuous. It seems to depend on the quark flavors and their masses. The extreme conditions 
in the early universe, i.e., high temperatures, high densities and out-of-thcrmal and out-of-chemical equilibrium, 
likely affect the properties of the partonic matter and control the dynamics of the phase transition. The equation 
of state deduced from lattice QCD calculations (and quasi particle model) plays a very essential role in present 
work. It sets the validity of the entire treatment. The high temperatures (energies), at which the strong coupling 
as nearly vanishes, defines the upper end of limitation. The lower one is characterized by the hadronic era. 
When applying Eckart theory, we find that the evolution of the Hubble parameter follows the same line defined 
by the standard cosmological model. The comparison with various types of matter shows that the comoving 
time behaves very smooth with H , although viscous hadron-QGP results in larger t than in viscous QGP. In 
both of them, t seems to be larger than in the collisionfree and ironviscous ideal matter. Israel-Stewart theory 
is assumed to solve the constrains of Eckart theory. Therefore, reliable results are to be expected. In order to 
make a qualitative estimation for the effect of viscosity, we compare the time evolution of the Hubble parameter 
in a viscous and nonviscous background matter. Apparently, we find that the viscosity drastically slows down 
the evolution. Should this result be confirmed, the whole picture about the evolution of early Universe has 
to be revised, accordingly. In order to make an estimation for this effect, the dynamics of phase transition(s). 
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interaction(s) and out-of-equilibrium and dissipative processes should be taken into account. The effect of the 
cosmological constant on the anisotropy and homogeneity and the cosmological density perturbations in the 
early Universe would play an essential role in characterizing the evolution of the cosmological parameters as 
well. 
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